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ABSTRACT 

We present the first radiation magnetohydrodynamics simulations of the atmosphere of white dwarf stars. We 
demonstrate that convective energy transfer is seriously impeded by magnetic helds when the plasma-/3 param¬ 
eter, the thermal to magnetic pressure ratio, becomes smaller than unity. The critical held strength that inhibits 
convection in the photosphere of white dwarfs is in the range B = 1-50 kG, which is much smaller than the typ¬ 
ical 1-1000 MG held strengths observed in magnetic white dwarfs, implying that these objects have radiative 
atmospheres. We have then employed evolutionary models to study the cooling process of high-held magnetic 
white dwarfs, where convection is entirely suppressed during the full evolution (B > 10 MG). We hnd that 
the inhibition of convection has no effect on cooling rates until the effective temperature (Tefi) reaches a value 
of around 5500 K. In this regime, the standard convective sequences start to deviate from the ones without 
convection owing to the convective coupling between the outer layers and the degenerate reservoir of thermal 
energy. Since no magnetic white dwarfs are currently known at the low temperatures where this coupling sig- 
nihcantly changes the evolution, effects of magnetism on cooling rates are not expected to be observed. This 
result contrasts with a recent suggestion that magnetic white dwarfs with Tj-ff < 10,000 K cool signihcantly 
slower than non-magnetic degenerates. 

Keywords: convection - magnetohydrodynamics (MHD) - stars: magnetic held - stars: evolution - stars: 
fundamental parameters - stars: interiors - white dwarfs 


1. INTRODUCTION 

Magnetic white dwarfs are stellar remnants featuring 
global magnetic structures with held strengths from 1 kG to 
1000 MG. They account for a signihcant part of the white 
dwarf population with an e stimated fraction of around 10% in 
volume complete samples (iLiebert et al.ll200^ ISchmidt et al.l 
120031 iKawka et aH 120071) . Most of these objects are high- 
held magnetic white dwarfs (HEMWD), with held strengths 
B > \ MG, and a distribution of mag netic held strengths 
that appears to pea k around ~ 20 MG (ISchmidt et al.llToO^ 
iKiilebi et al.l 1200^ . HEMWDs show obvious Zeeman line 
splitting in spectroscopic observations. It is currently diffi¬ 
cult to understand these data owing to the lack of an appro¬ 
priate theory of Stark broadening in the presence of a back¬ 
groun d magnetic held in an arbitrary direction (iMain et al.l 
Il998h . In particular, the standard spectroscopic technique 
emplo yed to derive atmosph eric parameters from the Balmer 
lines (IBergeron et al.l Il992h can not b e used to constrai n 
the mass and cooling age of HEMWDs (IKiilebi et al.ll200^ . 
However, there is a growing sample of HEMWDs in common 
proper motion pairs or with known trigonometric parallaxes, 
allowing to de rive masses. This sample shows a mea n mass 
of ^ 0.80 Mq (iFerrario et al.l2015l : lBriggs et al.l20r5h . which 
is signihcantly higher than the me an mass of non-magn etic 
white dwarfs 0.60 Mq, see, e.g.. lKleinman et ^120131) . 

Numerous recent studies have provided scenarios for the 
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origin of HEMWDs, accounting for their mass, velocity, and 
magnetic held strength distributions. The lack of a signih¬ 
cant trend in the number of HEMWDs as a function of at - 
mospheric composition and cooling age (IKiilebi et al.ll2009l) . 
as well the presence of held st rengths too large to be pro¬ 
duced by a convective dynamo (iDufour et al.ll200^ . suggest 
that magnetic helds are remnants of the white dwarf progen¬ 
itors. (Current scenarios tend to be grouped in three cate¬ 
gories, suggesting that HEMWDs are: rem nants of intermedi¬ 
ate mass stars with conserved foss il helds (I Angel et al.ll 19811 
Wickramasinghe & Ferrarioll2005h : the outcome of mergers 
of either two white dwarfs or a white dwarf and the core 
of a giant star (iGarcfa-Berro et al.ll2012UKulebi et al.ll2013al 

IWickramasinghe et al.ll2014l: iBriggs et al.ll2015l) : or products 
of the amplihcation of a seed held b y a convective dynamo in 
the core of the evolved progenito rs (iRuderman & Sutherland! 
119731 iKissin & Thompsoni 120151) . The origin of magnetic 
white dwarfs remains elusive since current observations do 
not allow to clearly differentiate between these evolution 
channels. Magnetic white dwarfs in clusters and common 
proper motion pairs (IKiilebi e t al.l 1201 01: iDobbie et al.l 120121: 
iKtilebi et al.ll2013bl: lo^ie et al.l UoTIb are consistent with 
single star evolution in some but not all cases. It is also 
difficult to explain the white dwarfs with the strongest mag¬ 
netic helds, as well as the absence of HEMWDs with late- 
type star companions, without invoking the merger sce nario 

(IWickramasinghe et akllToTil: IKissin & TTiompsonl2015l) . 

On the other hand, weaker magnetic helds (B < 1 MG) are 
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also found in white dwarfs, although they are difficult to de¬ 
tect systematically owing to the lack of obvious Zeeman split¬ 
ting i n high-resolution spectra for B < 20 kG (iJordan et alJ 
120071) . A few small spectropolarimetric surveys, sensitive 
to field strengths as small as kG, have put the fraction 
of kG-range white dwarfs at 3-30%, with a roughly con- 
stant 1-10% incidence per decade of magnetic field strength 
(iJordan et al.ll200^ iLandstreet et al.ll201^ iKawka & VennesI 

I 2 OI 2 I) . The uncertain ratio reflects the small number statistics 

in the current samples. These magnetic white dwarfs appear 
to have average mass es and are though t to origin from sin¬ 
gle stellar evolution (I.Tordan et alJl200'^ . with the magnetic 
fields possibly generate d through a dynamo process in the 
white dwarf progenitor (IWickramasinghe et al.ll201^ . Fur¬ 
thermore, large observed values for (Bz)/(|B|) are a strong 
indicator that the fields have a global organized structure, un¬ 
like t he complex magnetic fields at the surface of Sun-like 
stars (ILandstreet et al.l2012^ . 

Gaia will provide precise parallaxes for more than 100,000 
white dwarfs, including all known mag netic white dwarfs 
(iTorres et alJ 120051; ICarrasco et al.l 120141) and spectroscopic 
follow-ups will identify even more magnetic objects. Gaia 
will establish the first homogeneous mass distribution and 
cooling sequence of magnetic remnants. Given the ubiq¬ 
uitous presence of magnetic white dwarfs in the high-mass 
regime, it is critical to understand these objects to recover the 
Galactic star form ation history and initia l mass function in the 
^ 3-8 Mq range (iTremblav et al.llT014ft . Magnetic remnants 
can also be used to constrain st ellar evolution at intermedi¬ 
ate masses (iKiile hi e t al.l 1201 3bh and study possible popula¬ 
tions of mergers (iBadenes & MaozI 1201 2t IWegg & PhinnevI 
I 2 OI 2 I) . It is therefore essential, at this stage, to build pre¬ 
cise model atmospheres and evolution sequences for these 
peculiar degenerate stars. It has been suggested for a long 
time that convection is complete l y inhibited in HFM WDs 
dWickramasinghe & Martini 1 19861 : IValvavin et al.l 120141) . al¬ 
though this has not yet been ver i fied w ith realistic simula¬ 
tions. Furthermore, iKenler et alJ (1201 3l) suggest that small 
undetected magnetic fields could impact the mass distribution 
o f cool convective white dwarfs. 

IValvavin et al.l (l2014ft have recently proposed that the in¬ 
hibition of convection in magnetic white dwarfs has a large 
impact on cooling rates, by increasing the cooling times by a 
factor of two to three. However, they arrived at this conclusion 
with a simple analytical argument and it needs to be verified 
with state-of-the-art evolution models. In this work, we per¬ 
form the first radiation magnetohydrodynamics (RMHD) sim¬ 
ulations of the atmosphere of magnetic pure-hydrogen (DA) 
white dwarfs (Section 2.1). We then consider the results 
of these simulations for the computation of new cooling se¬ 
quences for magnetic white dwarfs using an established evo¬ 
lution code (Section 2.2). We discuss the implications of these 
results in Section 3 and conclude in Section 4. 

2. WHITE DWARF MODELS 
2.1. Magnetohydrodynamics Simulations 

We have computed RMH D simulations with the CO^BOLD 
code (iFrevtag et al.l 120121) for pure-hydrogen DA white 
dwarfs. We rely on a representative set of atmospheric pa¬ 
rameters, Teff ^ 10,000 K and a surface gravity of logg = 8.0, 
and our simulations are detailed in Table 1. The setup of the 
simulations is v ery similar to that of the non-magnetic mod¬ 
els presented in ITremblav et al.l (l2013allbl) . In the tempera¬ 


ture regime considered here, the convection zone is signifi¬ 
cantly deeper than the atmospheric layers, and we use a bot¬ 
tom boundary (at Rosseland optical depth tr ^ 10^) open to 
convective flows and radiation, where a zero net mass flux 
is ensured. We fix the entropy of the ascending material at 
2.0819 erg g“' K“' for all simulations, corresponding to the 
value in the non -magnetic simulation at = 10,025 K and 
logg = 8.0 from ITremblav et aP (l2013bl) . The lateral bound¬ 
aries are periodic, and the top boundary is open to material 
flows and radiation. We rely on the same opacities, equa¬ 
tion of state, and grid resolution (150 x 150 x 150) as the pre¬ 
vious non-magnetic simulations. The frequency dependent, 
i.e. non-gray, radiative transfer is solved along long charac¬ 
teristics employing Feautrier’s method. Opacities are grouped 


into bins (for details on the opacity binning approach see, e.g.. 

lNordlundll982l:lLudwig et al.lll994 

IVogler et al.ll2004l) using 

the 8-bin scheme oflTremblay et al 

(I2013bh. 


The main difference compared to earlier computations is 
that we have imposed, at the start of the simulations, a verti¬ 
cally oriented magnetic field (towards the exterior of the star) 
with an amplitude of 0.5 kG and 5 kG for our two magnetic 
simulations, respectively. The magnetic boundary conditions 
are imposed independently to the hydrodynamics conditions. 
We require that magnetic field lines remain vertical at both 
the top and bottom layers, while lateral boundaries are peri¬ 
odic. We further require that the magnetic flux is constant at 
the bottom, mimicking the effect of a global fossil field an¬ 
chored in the deep degenerate core. We note that our RMHD 
simulations do not assume hydrostatic equilibrium and auto¬ 
matically take into account the turbulent pressure, magnetic 
tension forces, and magnetic pressure. 

The MHD mod ule of CO^BOLD (see Section 3.7 of 
IFrevtag et alJl2012h provides several numerical methods for 
solving the MHD equations, which are quite different from 
the ones employed for pure hydrodynamics. In particular, 
the m ethod used here relies on the HLL solver (iHarten et alJ 
119831) . which is more stable but with increased dissipation 
compared to the Roe solver used for the published grid of 
DA white dwarfs. The MHD module also enforces the 
divergence-free condition V ■ B = 0 bas ed on a constrained 
transport scheme (see, e.g., lT6thll2000h . In order to study 
the effect of magnetic fields on the atmospheric stratification, 
we have computed a third model with the same MHD solver 
but no magnetic field. We computed all simulations for five 
seconds in stellar time, which is several times the convective 
turnover timescale. 

Figure [T] presents snapshots of the emergent intensity for 
our three relaxed simulations. From an average over 125 
snapshots, we also display at the top of the panels the 
values (derived from the emergent flux) and the relative in¬ 
tensity contrast. We observe that magnetic fields have a sig¬ 
nificant impact on the emergent intensity. For Bz = 0.5 kG, 
diverging upflows concentrate magnetic flux in downflows, 
much lik e what is observed in the so-called quiet regions on 
the Sun dNordlund et alj |2009|) , which are characterized by 
a rather weak average magnetic flux. Small magnetic flux 
concentrations form and appear as bright intergranular points 
since they act as radiative leak due to their reduced mass den¬ 
sity. Table 1 demonstrates that the root-mean-square verti¬ 
cal magnetic field in the photosphere is significantly larger 
than the average magnetic field owing to these flux concentra¬ 
tions. For a field strength of Bz = 5 kG, convection is already 
largely inhibited, and occurs as narrow and bright plumes 
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Table 1 

RMHD Models of White Dwarfs 


Tiff 

(K) 

logg 
(cm s“^) 

X 

(km) 

z 

(km) 

z(tR = l)-Zbot 
(km) 

Bz 

(kG) 

5z,rms('Tl — 1) 

(kG) 

/ (7) 
(%) 

Mach(TR = 1) 

10024 

8.0 

2.11 

0.83 

0.38 

0 

0 

14.69 

0.45 

10037 

8.0 

2.11 

0.83 

0.38 

0.5 

1.28 

14.13 

0.38 

9147 

8.0 

2.11 

0.63 

0.30 

5.0 

5.38 

21.86 

0.25 


Note. — All quantities were averaged over 250 snapshots and over constant geometrical depth when appro¬ 
priate. Tgff is derived from the temporal and spatial average of the emergent flux. Bz is the horizontally averaged 
magnetic field, which is constant at all times and all depths from the requirement of magnetic flux conservation. 
Sz,rms is the rms vertical magnetic field at th e geometrical depth that corresponds to (rR)x,y = 1- <5/nns/(f) is 
the relative intensity contrast (see Eq. (73) of IFrevtag et al.ll2Q12D . 


very similarly to Sun spots wh ere B > 2.5 kG (IWeiss et al.l 
Il996tl^hussler & Vogleiti2006h . This is not a surprising re¬ 
sult since the thermal pressure in the photosphere of the simu¬ 
lated white dwarfs is only slightly larger than that in the Sun, 
and a similar magnetic pressure is necessary to inhibit con- 
yectiye flows. Studies of the impact of magnetic helds on 
surface conyection i n the Sun and Sun-like stars by numerous 
RMHD simulation s (iRemp el et al.ll 2009t IChe ung et a l]l2010t 

iFreytag et al.l2012l:lBeeck et al.ll2013l:ISteiner et al.l2014b can 

also be used to learn about the same process in white dwarfs, 
eyen though the origin and large scale structure of magnetic 
helds are yery different. 

Figure |2] presents the temperature prohles of our simula¬ 
tions, drawn from the ayerage of {T'^) oyer surfaces of con¬ 
stant tr for 12 snapshots. For the 0.5 kG simulation, we ob- 
serye that the magnetic held only has an impact on the up¬ 
per photosphere (tr < 10“^), where the temperature gradient 
is shallower. The importance of the feedback effect of mag¬ 
netic helds on the stellar structure can be estimated from the 
plasma-/3 parameter 


the thermal to magnetic pressure ratio, where P is the thermal 
pressure, and B the ayerage magnetic held strength. Since 
the thermal pressure is rapidly decreasing with height while 
the magnetic pressure is roughly constant, magnetic feed¬ 
back effects increase with height. There are two main rea¬ 
sons for the shallower temperature gradient in the uppermost 
stable layers of magnetic white dwarfs. First of all, mag¬ 
netic held lines restrain conyectiye hows, hence the conyec- 
tiye oyershoot that usu ally cools the upper layers is weaker 
(iTremblay et alJ2013ah . This is a purely dynamical effect that 
will not be obseryed in a ID magnetic model with local con¬ 
yection. Furthermore, the radiatiye heating, magnetic dissipa¬ 
tion, and magnetic pressure all contribute to increase the ther¬ 
mal pressure scale height compared to the non-magnetic case, 
which implies a shallower temperature gradient as a function 
of geometrical depth. In general, the consequence is also a 
shallower temperature gradient as a function of tr. 

For a held strength of 5 kG, the oyerall atmospheric stratih- 
cation is signihcantly impacted by the presence of a magnetic 
held. Conyectiye energy transfer is impeded in the photo¬ 
sphere and Figure [3 demonstrates that the conyectiye hux at 
Tr = 1 is reduced by a factor of two compared to the non¬ 
magnetic model. The smaller conyectiye energy transfer im¬ 
plies that the stratihcation in the conyectiyely unstable regions 
must adjust to a steeper temperature gradient to transport the 
same amount of total Hux. The temperature gradient in the 


line-forming regions becomes yery close to the radiatiye gra¬ 
dient, as demonstrated in Figure|2]from the comparison with a 
ID structure where conyection was artihcially suppressed. On 
the other hand, in the deeper layers where the thermal energy 
is larger, conyection is still significant for this held strength. 
Neyertheless, the steeper temperature gradient in the upper 
conyectiyely unstable layers (tr > 0.1), caused by the inhibi¬ 
tion of conyection, decreases the Teff yalue by 880 K for the 
same conditions at the bottom. Full eyolutionary calculations 
are necessary to link the magnetic atmospheres to the stellar 
interior, and this result does not imply that magnetic white 
dwarfs haye smaller luminosities for the same core tempera¬ 
ture (see Section 2.2). For our models at Teff ~ 10,000 K and 
logg = 8.0, /3 = 1 for B ^ 5.7 kG at the photosphere (tr = 
1). This critical held strength is yery close to the obseryed 
transition between a conyectiye and an almost fully radiatiye 
temperature gradient in the RMHD simulations. Our results 
support the suggestion that when the plasma-/3 is smaller than 
unity, i.e. when the magnetic pressure dominates oyer the 
thermal energy, the white dwarf atmospheric stratihcation ad¬ 
justs to a radiatiye gradient since conyectiye energy transfer is 
signihcantly hampered. 

In those cases where the plasma-/3 parameter is smaller than 
unity, the atmosphere is not expected to become static or ho¬ 
mogeneous since the stratihcation is still conyectiyely unsta¬ 
ble, albeit unable to create energetically efficient conyectiye 
hows. In particular, the relatiye intensity contrast for the B = 
5 kG simulation is still 21.9%, an eyen larger yalue than for 
the non-magnetic simulation. While conyection is restricted 
to narrow and inefficient plumes, the temperature contrast and 
yelocities in those structures are still large. It is currently 
unclear how these huctuations would decrease as the mag¬ 
netic held strength is further increased. It is a serious techni¬ 
cal challenge to compute RMHD simulations with larger held 
strengths since the time steps are dictated by the Alfyen speed 
B j where p is the density. For instance, the simulation 

at 5 kG is already of the order of 10 times slower than the non¬ 
magnetic simulation. Finally, the magnetic held tends to form 
localized hux concentrations in the intergranular lanes, and 
the spatial resolution of our RMHD simulations likely needs 
to be improyed in order to properly characterize the intensity 
contrast and small-scale huctuations. 

We h aye employed a stand ard grid of ID model atmo¬ 
spheres (ITremblay et al.ll201 Ih to compute the critical mag¬ 
netic held strength, dehned by /? = 1, aboye which conyection 
is signihcantly suppressed in the photosphere (tr = 1). Fig¬ 
ure U shows that the critical held is always below ^50 kCj. 
Known magnetic white dwarfs haye held strengths typically 
much larger than these yalues, and our results suggest that 
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Figure 1. Bolometric intensity emerging from the xy plane at the top of computational domain for CO^BOLD 3D simulations computed with the MHD solver. 
All simulations have a constant surface gravity of logg = 8.0, a pure-hydrogen composition, and rely on the same entropy value for the inflowing material through 
the open bottom boundary. The simulations shown in the middle and right panels have at the bottom an imposed average vertical magnetic field of 0.5 kG and 
5 kG, respectively. The rms intensity contrast with respect to the mean intensity and Tgff values ai‘e also shown above the panels. The length of the bai‘ in the top 
right is 10 times the pressure scale height at tr = 2/3. 
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Figure 2. Temperature structure as function of optical depth (logarithmic 
scale) for the non-magnetic (black), 0.5 kG (blue), and 5 kG (red) 3D RMHD 
simulations. The temperature was determined from an average of {T‘^) over 
surfaces of constant optical depth. We also show a purely radiative ID model 
atmosphere (green), where convection was aifificially suppressed, at the same 
Teff value as the 5 kG model. 

convection is suppressed at the surface of HFMWDs. Fur¬ 
thermore, while we have only performed simulations with 
a vertically oriented magnetic held, it is generally thought 
that the damping of convection is even stronger for horizon¬ 
tally oriented helds since the Lorentz force will act against 
vertical Hows. In other words, convection is expected to 
be globally inhibited above a certain magnetic Held strength 
(IValvavin et al.ll20i4b . 

The rapid increase of ^ as a function of depth implies that 
when convection is suppressed at the surface, it could still be 
fully developed in deeper layers as demonstrated by our 5 kG 
simulation. Once /3 = 1 at the base of the convection zone, the 
entire convection zone is likely to be significantly disrupted. 
Figure 0] shows this critical field strength (dash ed lines) as 
predicted by ID envelopes (iFontaine et al.ll2()()Th . In the in¬ 
termediate regime between the suppression of convection at 
the surface and in the full convection zone, one should use 
radiative atmospheres but consider the possibility of an inter¬ 
nal convection zone. However, new cooling sequences with 
partial convective inhibition would need to be computed to 
determine the size and structure of these internal convection 


Figure 3. Ratio of the convective to total energy flux as a function of optical 
depth (logarithmic scale) for the non-magnetic (black), 0.5 kG (blue), and 
5 kG (red) 3D RMHD simulations. The (3D) convec tive flux is the sum of 
the enthalpy and kinetic energy fluxes (see Eq. 5 of ITremblav et al.ll2015l) 
averaged over constant geometrical depth. 


zones. These calculations are outside the scope of this work 
because a realistic magnetic field geometry would be required 
to properly model individual white dwarfs. Furthermore, it is 
difficult to extrapolate our RMHD results for the atmosphere, 
where convective velocities are close to the sound speed, to 
deeper convective layers where the flows have a kinetic en¬ 
ergy density that becomes far smaller than the thermal energy 
density. Once the magnetic field becomes larger than the ki¬ 
netic equipartition field strength 

</87r=i(pvLv), (2) 

where Vconv is the local convective velocity, different modes of 
convection with smaller physical scales may set in. Figure |5] 
demonstrates that the kinetic equipartition field strength is in 
the kG-range throughout the convection zone for a represen¬ 
tative 0.6 Mq white dwarf. It suggests that convection could 
be disrupted for magnetic field strengths smaller than those 
defined by the conservative P = I estimate of Figure]?] for the 
bottom of the convection zone. 
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Figure 4. Magnetic field strength that corresponds to plasma-/3 = 1 in the 
photosphere (tr = 1, solid lines) and base of the convection zone (dashed 
lines) as a function of Teff. Sequences are color-coded for logg = 7.0 (red), 
8.0 (blue), and 9.0 (black), from bottom to top. Plasma-/3 = 1 estimates when 
convective energy transfer is suppressed by the magnetic field. Photospheric 
values are derived from ID model atmospheres wi th a mixing-length param¬ 
eterization of ML2/« = 0.8 iTremblav et alj|701 ID while values for the base 
of the convection zone are derived from stand ard envelope models using a 
slightly more efficient ML2 /q: = 1.0 convection iFontaine et alj|200Th . 



Figure 5. Contours of kinetic equipartition magnetic field strength (loga¬ 
rithmic values in Gauss identified on the panel, see Eq. (2)) as a function of 
fractional mass log^ = log(l -M(r)/M*) integrated from the surface and Teff. 
We rely on the evolution sequence at 0.6 Mq with ML2/a = 1.0 convection 
(see Section 2.2). We also show the position of three atmospheric layers (tr 
= 0.1, 1.0, 10.0, from top to bottom in solid green lines) and the degeneracy 
boundary {rj = 0, solid red line). 

2.2. Evolutionary Models 

It has been known for a long time that superficial convection 
has no influence whatsoever on the cooling time until the base 
of the convection zone reaches into the degenerate reservoir 
of thermal energy and couples, for the firs t time in the cooling 
process, the surface with that reservoir (iTassoul et al.lll99u 
iFontaine et al.ll2()()l1) . The convective coupling occurs at Ti-ff 
values lower than 6000 K in white dwarfs, hence the suppres¬ 
sion of convection is not expected to impact cooling rates for 
war mer remnants. This arg ument contradicts the suggestion 
from IValvavin et aT] (120141 see Fig. 3a) that the suppression 
of convection changes the cooling rates and explains the ob¬ 
served temperature distribution of magnetic white dwarfs, for 


which their coolest bin is at T^ff = 6000 K. To demonstrate it 
quantitatively, this section presents new evolution sequences 
that we have comp uted with our state-of-the-a rt white dwarf 
evolutionary code (IFontaine et al.ll200Tl 12013h . To fully ap¬ 
preciate the results, we also review the important properties 
of white dwarf cooling in Section 2.3. 

We computed a standard 0.6 Mq sequence with a C/O core, 
a helium envelope containing 10“^ of the total mass, and a hy¬ 
drogen outer layer containing 10“^ of the total mass. In partic¬ 
ular, it takes into account superficial convection as it develops 
with time relying on t he so-called ML2 /q; =1.0 version of the 
mixin g-length theory (iBohm & Cassinellil 197 It ITassoul et'SI 
119901) . We have computed an additional sequence where con¬ 
vection is totally suppressed, thus mimicking the maximum 
possible effect of magnetic inhibition, e.g. for held strengths 
of 10 MG or larger according to Figure |4l Both sequences 
are presented in Figure |6] (left panel), where the solid curves 
refer to the normal sequence, while the dotted curves refer 
to the “magnetic” sequence. The location of convective cou¬ 
pling is indicated by the first dashed vertical segment from the 
left. This corresponds specifically to the model with the base 
of its convection zone first entering the degenerate thermal 
reservoir from above (the upper boundary of that reservoir is 
defined by a local value of the electron degeneracy parame¬ 
ter of p = 0, where ykT is the chemical potential of the free 
electrons). When convective coupling occurs, Teff = 5527 K 
and the cooling age is 3.13 Gyr. Above Teff = 5527 K, there is 
no significant difference whatsoever between the behaviors of 
the two sequences, meaning that magnetic inhibition of su¬ 
perficial convection does nothing to the cooling process in 
this hotter phase. We have also computed sequences at 1.0 
Mq which are likely more representative of the HFMWDs. 
Figure |6] (right) demonstrates that the behavior is similar to 
the lower mass case, and convective coupling takes place at a 
only slightly higher temperature. 

Our evolutionary sequences demonstrate that the cooling 
rates, hence the relation between core and surface tempera¬ 
ture, must remain the same for magnetic and non-magnetic 
white dwarfs. We now try to reconcile this fact with the pre¬ 
diction from our RMHD simulations indicating that the in¬ 
hibition of convection by a magnetic held creates a steeper 
(radiative) temperature gradient in the outer convectively un¬ 
stable layers. Figure [T] presents the temperature profile of a 
model at ^ 6200 K from the standard evolution sequence 
at 0.6 Mq, along with the case where convection was sup¬ 
pressed for the entire cooling process. It confirms that even 
though there is a much steeper gradient at the surface of mag¬ 
netic white dwarfs, this is not the case for all internal layers, 
and the non-magnetic relation between core and (average) sur¬ 
face temperature holds. Interestingly, Figure [8] demonstrates 
that for the magnetic case, the steep radiative gradient in the 
outer layers is associated with a very sharp opacity peak as a 
function of fractional mass. It is unclear if such opacity peak 
could generate pulsations in magnetic white dwarfs, which we 
discuss in Section 3.4. 

2.3. The Cooling Process in White Dwarfs 

We have designed Figure |9] to review the cooling process 
in white dwarfs. The cooling time depends on the amount 
of thermal energy contained in the star and the rate with 
which this energy is transferred from the thermal reservoir 
to the surface. The available thermal energy at a given 
epoch is given by the integral shown on the y-axis of Fig- 
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Figure 6. (Left:) cooling sequences in terms of cooling time (black, left axis) and central temperature (red, right axis) as a function of decreasing Teff for a 
0.6 Mq da white dwarf (solid lines). We have assumed thick helium and hydrogen layers with fractional masses of 10'^ and 10^, respectively. We have also 
computed a sequence where convection was artificially suppressed mimicking the effect of a strong magnetic field (dotted lines). Convection has no effect on 
the cooling until there is a convective coupling with the degenerate core at the position illustrated by the first dashed vertical segment from the left. The location 
where the cooling times become larger for the magnetic sequence is indicated by the second dashed vertical segment from the left. The Teff values for both 
transitions are shown on the panel. (Right:) same as left panel but for a 1.0 Mq DA white dwarf. 




Figure 7. Temperature structure as a function of fractional mass (both log¬ 
arithmic values) for a DA white dwarf at Teff ~ 6200 K and 0.6 Mq. The 
solid sequence (Tgff = 6243 K) relies on ID convection (ML2 /q: = 1.0) while 
the dotted sequence (Tgff = 6205 K) had convection suppressed in the entire 
cooling process. 

ure |9] Here, we show the running integral (black dotted 
curve), from the center to the surface, for three models be¬ 
longing to the standard (convective) evolutionary sequence at 
0.6 Mq discussed above. The x-axis shows the fractional 
mass logq' = log(l -M(r)/M*) integrated from the surface. 
The upper boundary of the reservoir of thermal energy is a 
concept that is a bit fuzzy, but it must correspond to a loca¬ 
tion on the flat part of each curve, i.e., to a layer above which 
there is practically no more contribution to the reservoir. Con¬ 
veniently, this boundary is usually dehned as the layer where 
the degeneracy parameter p = 0. In Figure |9] the layer ry = 0 
corresponds, for each model, to the location of the sharp cut¬ 
off on the left of the blue spike. With cooling, the boundary 
77 = 0 moves up toward the surface because the star becomes 


Figure 8. Similai' to Figure [T] but for the Rosseland mean opacity (ft) as a 
function of fractional mass (both logarithmic values). The solid sequence 
(Teff = 6243 K) relies on ID convection (ML2/a = 1.0) while the dotted se¬ 
quence (Teff = 6205 K) had convection suppressed in the entire cooling pro¬ 
cess. 

globally increasingly more degenerate. Moreover, we have il¬ 
lustrated in red the profile of the ratio of the convective flux to 
total flux, Fconv/Ttot. It should be understood that convective 
coupling arises when the base of the convection zone reaches 
the boundary 77 = 0, which is imminent but has not yet oc¬ 
curred in the coolest model (Teff = 5585 K) shown in the plot. 
In this particular evolutionary sequence, convective coupling 
occurs when the star has cooled down to the somewhat lower 
value of Teff = 5527 K. 

In a cooling white dwarf, the degenerate core and reservoir 
of thermal energy is relatively well insulated by a nondegener¬ 
ate envelope whose global opacity regulates the rate of energy 
loss. To illustrate this opacity barrier, and in particular the role 
of the insulating layers between the base of the outer convec- 
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0. -4. -8. -12. -16. -20. 

log q = log (1-M(r)/M.) 

Figure 9. Logarithmic value of the available thermal energy integrated from 
the center (left to right, black dotted curves) at three given epochs of the 
standard 0.6 Mq cooling sequence (Teff = 14011 K, 9752 K, and 5585 K, 
identified on the panel) within a certain fractional mass log^. The uppermost 
degenerate layer {rj = 0) corresponds, for each model, to the location of the 
sharp cutoff on the left of the blue spike. The blue spikes correspond to the 
running integral of the optical depth, from the base of the convection zone on 
the right to the layer where 77 = 0 on the left, with a scale of tr = 5 x 10^. The 
ratio of the convective to total flux is illustrated by the red profiles (ML2 /q; 
= 1.0) and the two black dots on each curve indicate, respectively, the depth 
where the magnetic pressure is equal to the gas pressure assuming a magnetic 
field of 10 MG (on the left) and 1 MG (on the right). We also indicate the 
location of three atmospheric layers along the x-axis, corresponding to tr = 
10, 1, and 0.1, from left to right (the short vertical dotted line segments). 

tion zone and the reservoir, we integrated the optical depth 
drR = -npdr between the base of the convection zone and the 
layer rj = 0. For each model considered, we plotted in Fig¬ 
ure |9] in blue the running integral of optical depth, from right 
to left together with a scale of tr = 5 x 10®. The blue spikes 
thus identify the layers that are of importance in the insulating 
process and in the role of regulator of the rate of energy trans¬ 
fer from the core to the surface. Even for the coolest model 
illustrated here, the opacity barrier is still enormous and the 
reservoir remains relatively well insulated. The convective 
coupling will occur in a somewhat cooler phase for which the 
base of the convection zone finally reaches the boundary ry = 0. 
From that point on in time, the reservoir becomes effectively 
coupled directly to the atmospheric layers via a convection 
zone whose efficiency reaches practically 100%. For the first 
time in the evolution of the star, the exact physical conditions 
characterizing the atmospheric layers will start playing a role 
in the cooling process. 

The layers where the blue optical depth curve is flat in Fig- 
ure|9]have a negligible contribution to the opacity barrier and, 
thus, cannot play any role in the cooling process. For ex¬ 
ample, for the two warmest models, all of the layers above 
log(l -M(r)/M*) ^ -9 have no impact on this process. For 
the coolest model, all of the layers above log(l -M(r)/M*) ^ 
-7 have no impact either; the insulating layer represented by 
the small blue spike is still extremely efficient at regulating 
by itself the outflow of energy. In this context, we have added 
in the figure two black dots on each curve which indicate, 
respectively, the depth where the magnetic pressure is equal 
to the gas pressure assuming a magnetic field of 10 MG (on 


the left) and 1 MG (on the right). These layers sit far above 
the opacity barrier, hence magnetic effects, namely magnetic 
pressure, may impact the actual stratification of these outer 
layers, but these layers play no role in the cooling process. 
They have negligible contribution to the energy reservoir and 
negligible contribution to the opacity barrier. 

A last view on convective coupling can be made from 
Figure [To] with the standard convective cooling sequence at 
0.6 Mq. The small dots represent the opacity contours, while 
the bold dots represent the convective layers. The opacity 
maximum is caused by hydrogen recombination. We also 
show the position of the degeneracy boundary (77 = 0 ) with 
a solid red curve. It is observed that when the degener¬ 
acy boundary crosses the convection zone, there is a radical 
change in the envelope stratification, and conductive transfer 
dominates for regions below the degeneracy boundary. Fig- 
ure| 6 ]also demonstrates that forboth the 0.6 and 1.0 M 0 cases, 
the cooling time of the normal convective sequence becomes 
larger than that of the magnetic sequence in the phase follow¬ 
ing the onset of convective coupling, while the central tem¬ 
perature immediately drops below th at of the magnetic mo del. 
This behavior has bee n explained bv lTassoul et al.l (Il990l) and 
iFontaine et al.l (1200 Ih . and it is perhaps best understood with 
the analogy of a warm oven. Convective coupling is like open¬ 
ing the door of the oven; there is initially an excess of heat 
coming out of the oven, while the inside temperature drops 
immediately. In a white dwarf undergoing convective cou¬ 
pling, the excess of thermal energy is translated into a delay 
in the cooling process and the cooling time increases accord¬ 
ingly. After this excess energy has been radiated away, con¬ 
vective coupling enters a second phase, and that is that of ac¬ 
celerated cooling because convection now couples for good 
the energy reservoir and the surface, and it transfers energy at 
a greater rate than radiation alone could do. It is thus only in 
this second phase of the process that the cooling time of the 
magnetic sequence becomes larger th an the cooling time o f 
the normal sequence, as suggested by IValvavin et al.l (120141) . 
In Figure | 6 | the vertical dashed line segments, marked Teff 
= 3340 and 3840 K for the 0.6 and 1.0 Mq models, respec¬ 
tively, indicate the very low Teff values below which this sec¬ 
ond phase can proceed. 

We conclude this section with a comparison to the cooling 
process in magnetized neutron stars, which is also regulated 
by a heat blanketing envelo pe between the atmo sphere and 
the stellar interior (see, e.g.. jPotekhin et al']l2005h . For these 
objects, thermal conduction is the dominant energy transfer 
mechanism in the degenerate electron gas within the insulat¬ 
ing layers, and it has been established that the suppression of 
thermal conduction in the direction transvers e to the magnetic 
field lines can influe nce the cooling rates (lHernauislill985l: 
iPotekhin et al.ll200^ . In a white dwarf, however, the insu¬ 
lating region is non-degenerate and thermal conduction only 
takes place in the stellar interior, where changes in the con¬ 
duction rates are unlikely to impact the cooling process. Av¬ 
erage magnetic fields are also much weaker in white dwarfs 
in comparison to magnetized neutron stars. 

2.4. Magnetic Effects on Structures 

Figure [TT] compares the gas and magnetic pressure for a 
characteristic structure at 0.6 Mq and T^ff ^ 9750 K. We as¬ 
sume a 10 MG field at the surface and a conservation of the 
magnetic flux AirBr^ in the interior. This is obviously a rough 
description of the actual magnetic geometry in the interior, 
which is poorly constrained by observations. Nevertheless, 
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Figure 10. Opacity contours (small dots, logarithmic values in cm^ g“' iden¬ 
tified on the panel) in the envelope of a 0.6 Mq white dwarf as a function of 
fractional mass (logarithmic scale) and T^ff. We rely on the standard con¬ 
vective sequence with ML2/a = 1.0, and the convection zone is illustrated 
by the bold dotted contours. We also show the position of three atmospheric 
layers (tr = 0.1, 1.0, 10.0, from top to bottom in solid green lines) and the 
degeneracy boundary {r} = 0, solid red line). 



Figure 11. Thermal pressure profile (solid line) of a DA white dwarf struc¬ 
ture at 0.6 Mq and Teff = 9752 K, the same model as on the middle panel of 
Figure|3 We also show the magnetic pressure (dotted line) for a field strength 
of 10 MG at the surface and assuming the conservation of the magnetic flux 
in the interior. 

it demonstrates that magnetic effects could only play a role 
in the outer layers and at the very center, although there is 
no evidence that magnetic field lines reach the central re¬ 
gion. For the illustrated model, a fractional mass depth of 
log(l -M(r)/M*) = -9 corresponds approximately to a frac¬ 
tional radius of log(l -r/R^) = -2.3. Thus, magnetic fields 
(at the 10 MG level) could at best only in the outermost 
0.5% of the radius have an influence on the structure of these 
representative white dwarf models. As a consequence, we 
conclude that current mass-radius relations for non-magnetic 
white dwarfs will hold for magnetic remnants as well. 

3. DISCUSSION 

The computation of RMHD simulations for DA white 
dwarfs confirms that convective energy transport is seriously 


impeded by magnetic held lines when the plasma-/3 parameter 
is smaller than unity. As a consequence, radiative ID model 
atmospheres can be employed for magnetic white dwarfs with 
B > 50 kG according to Figure |4] The main shortcoming in 
the modeling of most known magnetic white dwarfs remains 
the spectral synthesis of th e Baimer lines accounting for bot h 
Stark and Zeeman effects ( jWickramasinghe & Martinill986h . 

3.1. Photometric Variability of Magnetic White Dwarfs 

It is difficult to explain from our results the large number of 
magnetic white dwarfs that show ph otometric variations of a 
few percent over their rotation period (iBrinkworth et al J2013l 
iLawrie et al.ll2013t IValvavin et al.ll^l4ft . We have demon- 
strated in Section 2.2 that the partial or total suppression of 
convection is unable to change the average surface tempera¬ 
ture until there is coupling between the convection zone and 
the degenerate core at low Teff values. As a consequence, we 
can not naturally explain global emergent intensity variations 
for the Teff values of known magnetic white dwarfs. Flowever, 
we note that if t he magnetic he l d is m oving at the surface, as 
hypothesized bv IValvavin et alJ (1201 Ih for WD 1953-011, the 
envelope would take some time to adjust to the new surface 
conditions. The Kelvin-Helmholtz timescale of the portion 
of the envelope including the entire convection zone is one 
estimate for this thermal relaxation time, which varies from 
about one second at 12,000 K to about 1000 yr at 6000 K. 
Since the cooling rates must remain constant according to our 
evolutionary models, the flux fluctuations created from this 
mechanism would average out over the full surface but not 
necessarily over the apparent stellar disk. 

We note that photometric variations are observed in 
hot magnetic white dwarfs where no convection is pre¬ 
dicted, hence it is already clear that convective effects 
are not involved in some cases. Previously supplied ex¬ 
planations for photometric variations remain valid, such 
as magneto-optical effects involving radiative transfer un¬ 
der different polarizations (Martin & Wickramasinghd ll979l : 
Wickramasinghe & Martini 1986nFerrario et al.ni997 ). Fi¬ 
nally, variations are also observed, although with a weaker 
amplitude, in apparently non-magnetic white dwarfs, where 
accretion hot spots or UV flux fluctuations and fluorescent 
optic al re-emission hav e been suggested as possible explana¬ 
tions (iMaoz et al.ll2015h . 

3.2. Cooling Age Distribution of Magnetic White Dwarfs 

Our results do not support the hypothesis that the observed 
distribution of HFMWDs as a function of Teff can be ex¬ 
plained by different cooling timescales between magnetic and 
non-magnetic white dwarfs. This does not imply that the 
number ratio of magnetic to non-magnetic remnants should 
be constant as a function of Teff. The cooling age distri¬ 
bution of HFMWDs could be different from the fact alone 
that they have a distinct mass distribution. A variation of 
the velocity distribution as a function of both mass and Teff 
(IWegg & Phinnevl2012h . a consequence of the different main- 
sequence lifetimes, could change the magnetic incidence as 
a function of Teff even for volume-complete samples. Fur¬ 
thermore, a distinction between magnetic and non-magnetic 
objects could be present if a significant fraction of magnetic 
white dwarfs originate from mergers, which presumably have 
a different cooling history compared to single remnants. Fi¬ 
nally, very cool DA white dwarfs have deep convection zones, 
and for Teff < 6000 K, they reach a regime where the convec¬ 
tive turnover timescale at the base of the convection zone is of 
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the order of a few hours, which is similar to the rotatio n pe¬ 
riods of magnetic white dwarfs (iBrinkworth et al.ll2013l) . The 
hypothesis of a aw convective dynamo becomes tantalizing, 
although this needs to be tested with dynamical models. How¬ 
ever, this dynamo is unlikely to generat e fields stronger than 
the kinetic equipart i tion field strength (Fontaine et al.l 119731 
iThomas et alJll9931 iDufour et al.l 12008 ). Figure |5] demon¬ 
strates that for our standard evolutionary sequence at 0.6 Mq, 
the equipartition field strength reaches a maximum value of 
Beq ^10 kG at the base of the convection zone, suggesting it 
is an unlikely scenario for the known magnetic white dwarfs. 

We have found no firm evidence in the literature for a vari¬ 
ation in the incidence of magnetic white dw arfs as a func¬ 
tion o f Teff, which differs from the claim of IValvavin et alJ 

(1201 4|) t hat the picture has n ow been settle d. _ On the 

contrary, iLiebert et alJ (l2003h . [Hollands et al.l (120151) . and 
iFerrario^ alJ (I2OI5I) suggest that variations still need to 
be confirmed owing to several observational biases and 
conflicting r e sults. Furthermore, iKlilebi et al.l (120091) and 
iKenler et al.l (1201 3h And no clear evidence of variations in 
the homogeneous SDSS sample, although mo st objects have 
Tj-ff > 7000 K, above the temperature where IValvavin et alJ 
(I2014^ observe a significant increa se. There is marginal evi- 
dence from the local 20 pc sample (iGiammichele et aO2012^ 
that the incidence of magnetic fields increases for Teff < 
6000 K. If we consider only DA white dwarfs as well as 
objects with a der i ved di stance under 20 pc in Table 2 of 
IGiammichele et all (120121) . we And a magnetic incidence of 
22± 11% (4 magnetic objects) for 5000 < Teff (K) < 6000, 
while the value is 10 ± 5% for warmer objects. We believe it 
is necessary to confirm this behavior with larger samples to 
fully understand the evolution of magnetic white dwarfs. 

3.3. Magnetic Fields in the White Dwarf Population 

Few magnetic white dwarfs have precise atmospheric pa¬ 
rameter determinations, and it is typical to exclude them from 
the samples employed t o derive the mean pro perties of held 
white dwarfs (see, e.g.. iTremblav et £01201 ll) . It is however 
difficult to detect magnetic objects with B < \ MG at low 
spectral resolution, hence it is therefore nearly impossible to 
define clean non-magnetic samples. 

We have shown that magnetic fields of a few kG can signif¬ 
icantly impact the thermal stratification in the upper layers of 
convective DA white dwarfs. Yet these fields are too weak 
to produce any significant Zeeman splitting, hence white 
dwarfs harboring s uch fields would not easily be detected. 
iKepler et alJ (120131) have suggested that undetected magnetic 
fields could explain the so-called high-l ogg problem observe d 
in the white dwarf mass distributions (iBergeron et alJ[T990l) . 
On the other hand, it was recently demonstrated that this 
problem is instead caused by inaccuracies in the I D mixing- 
length convection mo^l (ITremblav et all l2013bl) . Further¬ 
more, [^^lereLiDdlQlll) suggest that held strengths increase 
for convective objects, which would be a manifestation of the 
amplification of magnetic fields by convection. However, all 
their observations have B > \ MG, which is too strong to be 
amplified by convection since the kinetic equipartition held 
strength is always much smaller than B = I MG as demon¬ 
strated in Figure |5] 

To understand the effects of a population of white dwarfs 
with small undetected magnetic fields, we have computed 
synthetic ID spectra at Teff = 10,000 K and logg = 8.0. The 
spectra are derived from both a standard convective model at¬ 
mosphere with ML2/a = 0.8 (ITremblav et alJl^l ll) . and a ra- 




AA (A) AA (A) 

Figure 12. Left: Predicted Balmer line profiles for a DA white dwarf at 
= 10,000 K and log^ = 8.0. The spectra were computed from a standard 
ID model atmosphere with ML2 /q; = 0.8 convection (red) and a radiative ID 
atmosphere where convection was artificially suppressed (black), represent¬ 
ing the effect of a B > 5 kG magnetic field. All profiles are normalized to 
a unity continuum and the transitions are identified on the panel. We have 
employed a convolution of 3 A (FWHM) to represent typical observations. 
Right: Similar to left panel but we show the standard convective ID model 
spectra (red, Teff = 10,910 K and logg = 8.08) that best fits the radiative 
model. It suggests that fitting this magnetic white dwarf with a proper ra¬ 
diative model would result in Alogg = -0.08 and ATeff = -910 K compared 
to the standard convective solution. 


diative atmosphere where convection was completely inhib¬ 
ited, mimicking the effect of a weak 5 <B (kG) <100 mag¬ 
netic held, i.e. the range where Zeeman splitting is negligible 
at low spectral resolution. The left panel of Figure [T2ldemon- 
strates that the predicted Balmer lines of the two models are 
significantly different, although when projecting the magnetic 
model on a grid of convective models on the right panel of 
Figure [12] the Balmer lines look much alike, albeit with an 
offset in the atmospheric parameters. It implies that it would 
be difficult to identify such a small magnetic held from spec¬ 
troscopy alone. This could have an impact on the observed 
mass distribution of cool convective white dwarfs, although 
the logg shift is moderate according to Figure[T2| and the in¬ 
cidence of magnetic white dwarfs in the ^10 kG-range is ex¬ 
pecte d to be small (iKawka & Vennesll2012t iLandstreet et alJ 

iMl . 

The situation is different when accounting for 3D effects. 
In that more realistic case, the magnetic held inhibits convec¬ 
tive overshoot so that the upper layers (tr < 10<) must be in 
radiative equilibrium (see Figure|2]i. In ID, these upper layers 
are never convective and always in radiative equilibrium. As a 
consequence, synthetic spectral line cores based on 3D simu¬ 
lations are significantly shallower in the magnetic case, while 
they do not change in ID. We refrain from a quantitative pre¬ 
diction at this stage since the 3D RMHD structures have been 
computed with different numerical parameters in comparison 
to the published 3D grid. Neverthel ess, it is a potent i al expla - 
nation for the problem observed in ITremblav et al.l (l20131j) . 
where the predicted 3D line cores are systematically too deep 
and suggest that the upper layers are too cool. Figure |2| illus¬ 
trates that a held strength of ^ 1 kG is sufficient to signifi¬ 
cantly heat the upper layers. 

It is unlikely that the commonly proposed evolution sce¬ 
narios for magnetic white dwarfs could systematically gen¬ 
erate ^ 1 kG magnetic fields, which would then impact 
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the observed line cores. A plausible alternative, however, 
is that a turbulent dynamo systematically generates weak 
magnetic fields in convective white dwar fs, a well dis¬ 
cussed scenario for quiet regions of the Su n (ICattaneoll 19991 : 
IVogler & Schussl^lioOTirMoll et aO201 Ih . It consists of the 
amplification of small seed magnetic fields by the electrically 
conducting turbulent convective flows at the surface. Such 
fields are likely to reach an equilibrium strength of a fraction 
of the kinetic equipartition energy, corresponding to 0.1-1 kG 
in the photosphere of convective DA white dwarfs according 
to Figure|5] The magnetic fields would have characteristic di¬ 
mensions of the convective eddies of at most a few hundred 
meters, hence it would be difficult to detect them, except from 
their systematic feedback effect on the atmospheric stratifica¬ 
tion. As a consequence, recent spectropolarimetric surveys 
provide no direct constraint on this scenario. We hope to 
compute turbulent dynamo RMHD simulations in the future 
to verify whether the magnetic fields reach a sufficient ampli¬ 
tude to solve the discrepancy between the predicted 3D line 
cores and observations. 

3.4. Pulsating White Dwarfs 

It is difficult to apply our results quantitatively to pulsat¬ 
ing white dwarfs. The base of the convection zone corre- 
spon ds to the driving region of the ZZ Ceti pulsations (see, 
e.g.. iFontaine & BrassardI |2008^ . hence the dashed lines in 
Figure |4] illustrate the critical field where convective energy 
transfer will be largely suppressed in these layers. Thus, mag¬ 
netic fields stronger than 1 MG will likely have a dramatic 
effect on the driving mechanism of the pulsations, although 
it is difficult to rule out pulsating instabilities at this stage 
since the stratification will still be unstable. Another aspect 
of the problem is that the inhibition of convection will create 
a strong temperature gradient and opacity peak in the convec- 
tively unstable upper layers (see Figure [S]), which could inde¬ 
pendently drive pulsations through a K-type mechanism. This 
process has already been sugge sted for pulsating an d strongly 
magnetic hot DQ white dwarfs (iDufour et al.ll200§) . 

It is difficult to predict the position of an instability strip 
for magnetic white dwarfs since it is likely to depend on the 
strength and geometry of the magnetic fields. Indeed, mag¬ 
netic pressure will impact the position of the opacity peak as 
function of the radius. We n ote that the Lorentz force affects 
nonradial pulsations as well (ISaioll201^ . requiring additional 
theoretical work to model pulsating magnetic white dwarfs. 
However, the Ohmic timescale in the outer layers is short, 
suggesting that the magnetic field could be relaxed to a force- 
free potential state. This further highlights the fact that one 
must rely on realistic magnetic field geometries to model pul¬ 
sations in magnetic white dwarfs. 

For DA atmospheres, MG-range fie lds are excluded for the 
56 bright ZZ Ceti white dwarfs in the lGianninas et af] (1201 Ih 
sample, suggesting that magnetic fields inhibit pulsations. On 
the ot her hand, none of the HFMWDs with known and 
logg (iBriggs et alJl2015l) are within the ZZ Ceti instability 
strip, an essential ingredient to conclude about the possibil¬ 
ity of HFMWD ZZ Ceti white dwarfs. 

4. CONCLUSION 

We have computed the first RMHD simulations of pure- 
hydrogen white dwarf atmospheres. We have demonstrated 
that convective energy transfer is largely suppressed in the at¬ 
mosphere of magnetic white dwarfs for field strengths larger 
than B ^ 50 kG, confirming quantitatively the previously 


widespread idea that HFMWDs have no surface convection. 
Stronger magnetic fields are necessary to fully suppress con¬ 
vection in the envelope, and we find that for B = 1-100 MG, 
depending on the atmospheric parameters, the full stratifica¬ 
tion becomes radiative. For intermediate field strengths, the 
suppression of convection in the upper layers will change the 
stellar structure in a complex way, and new calculations with 
partial convective inhibition and realistic magnetic field con¬ 
figurations must be performed to better understand these ob¬ 
jects. 

We have presented new evolutionary calculations for DA 
white dwarfs where convection was fully suppressed, e.g. 
mimicking the effect of a B > 10 MG field. We find that 
the suppression of convection has no impact on the cooling 
rates until there is a convective coupling between the convec¬ 
tion zone and the degenerate core in the standard sequence 
at Teff ~ 5500 K. The currently known magnetic remnants, 
which are almost all above this temperature, are thus cooling 
like non-magnetic white dwarfs. Our results also suggest that 
the effect of magnetic pressure on the mass-radius relation is 
at most of the order of 1%. Finally, we conclude that the pho¬ 
tometric variations observed in a large fraction of magnetic 
white dwarfs remain largely unexplained. 
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